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ABSTRACT 



Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are com- 
piled from the literature. By empirical fitting we identify eight possible emission 
components and summarize the results in a "synthetic" lightcurve. Both optical 
flare and early shallow-decay components are likely related to long-term central 
p ■ engine activities. We focus on their statistical properties in this paper. Twenty- 

Q ■ four optical flares are obtained from 19 GRBs. The isotropic R-band energy 

is smaller than 1% of -E^^so- The relation between isotropic luminosities of the 
c3 ■ flares and gamma-rays follows L^^-^^o ^ -^7,^0''°^^ • Later flares tend to be wider 

and dimmer, i.e., w^ ~ tp/2 and -^^j^^ oc [tp/(l + ^)]-i-i5±o.i5_ rpj^^ detection 
probability of the optical flares is much smaller than that of X-ray flares. An 



> 

CN ■ optical shallow decay segment is observed in 39 GRBs. The relation between the 



break time and break luminosity is a power-law, with an index of —0.78 ± 0.08, 
CN ■ similar to that derived from X-ray flares. The X-ray and optical breaks are usu- 

^ ■ ally chromatic, but a tentative correlation is found. We suggest that similar to 

^ ■ the prompt optical emission that tracks 7-rays, the optical flares are also related 

to the erratic behavior of the central engine. The shallow decay component is 
likely related to a long-lasting spinning-down central engine or piling up of flare 
rS • materials onto the blastwave. Mixing of different emission components may be 

c^ ■ the reason of the diverse chromatic afterglow behaviors. 
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1. Introduction 

Gamma-Ray Bursts (GRBs) and their broadband afterglows are the most luminous 
phenomena in the Universe. According to the standard model, the broad-band afterglow is 
from the external shock as the fireball is decelerated by the ambient medium (Meszaros & 
Rees 1997; Sari et al. 1998). The prompt gamma-ray emission, on the other hand, is due 
to some internal dissipation processes within the relativistic ejecta, either due to internal 
shocks (Rees & Meszaros 1994; Kobayashi et al. 1997; Daigne & Mochkovitch 1998) or 
internal magnetic energy dissipation processes (e.g. Usov 1992; Thompson 1994; Drenkhahn 
& Spruit 2004; Giannios & Spruit 2006; Zhang & Pe'er 2009; Zhang & Yan 2011), either near 
or far above the photosphere. In the pie-Swift era, afterglow observations were mostly made 
in the optical bands. The data are well explained by the external shock model (e.g., Meszaros 
& Rees 1997; Sari et al. 1998; Panaitescu et al. 1998; Panaitescu & Kumar 2001; Huang 
et al. 2000; see Zhang & Meszaros 2004 for review). The successful launch and operation 
of the Swift mission (Gehrels et al. 2004) have significantly improved our understanding 
on the physical origin of GRBs. In particular, early X-ray afterglow observations revealed 
erratic flares and early plateaus that are difficult to interpret within the standard theoretical 
framework (Zhang et al. 2006, Nousek et al. 2006). The flares are believed to be produced 
by late activity of the GRB central engine (Burrows et al. 2005; Fan & Wei 2005; Zhang 
et al. 2006; Dai et al. 2006; Proga & Zhang. 2006; Perna et al. 2006), and the shallow 
decay segment likely signals a long-lasting wind powered by the GRB central engine after the 
prompt gamma-ray phase (Zhang et al. 2006). These features indicate that the GRB central 
engine does not die quickly. One is obliged to accept a more complicated afterglow picture, 
namely, the observed afterglow emission is a superposition of the traditional external shock 
afterglow and an afterglow related to the late central engine activity (Zhang 2011). 

The prompt localization of GRBs with the X-Ray Telescope (XRT) on board Swift 
has significantly increased the number of GRBs with optical afterglow detection and red- 
shift measurement. The simultaneous observations with XRT, UVOT, and ground-based 
optical telescopes in a much wider time window in the afterglow phase have revolutionized 
our knowledge about the GRB afterglow, and raised some critical problems with the con- 
ventional models (e.g., Zhang 2007, 2011; Liang 2010). For example, the multi-wavelength 
observations sometimes revealed chromatic behaviors between the optical and X-ray bands 
(e.g. Panaitescu et al. 2006; Liang et al. 2007), suggesting multiple emission components 
in the afterglow that dominate different energy bands at different epochs. The existence 
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of engine-powered early X-ray emission also obscures the clear separation of long and short 
GRB, causing confusion in GRB classification (e.g. Zhang et al. 2007, 2009; Lii et al. 2010 
for detailed discussion). In order to unveil the underlying physics, it is essential to decompose 
the lightcurve into different components that have distinct physical meanings. 

There are two approaches to decompose the lightcurves into different components: one 
through theoretical modeling and the other through empirical fitting. Theoretical modeling 
gives insights into the physical properties of the emission, including radiation mechanisms, 
micro-physical parameters, properties of the surrounding medium, etc. Intense modeling of 
optical afterglow data has been carried out in the pre- Swift era (e.g., Panaitescu et al. 1998; 
Panaitescu & Kumar 2001, 2002; Huang et al. 2000; Wu et al. 2005). This becomes increas- 
ingly difficult in the Swift era. First, there are many more bursts with high quality data to be 
modeled, and performing such modeling would be time consuming. More importantly, data 
suggest that it is essentially impossible to interpret all the data with the standard external 
shock model. The empirical fitting approach, on the other hand, makes use some empirical 
functions to fit the data, and is suitable to handle a great amount of data. A morphological 
study can catch insights of various emission components, which can ease statistical analyses 
of a large sample of data. This has been done by some authors (e.g., Liang & Zhang 2006; 
Panaitescu & Vestrand 2008, 2011; Kann et al. 2010; 2011). Some interesting features have 
been revealed. For example, Liang & Zhang (2006), Nardini et al. (2006), and Kann et al. 
(2006) found two universal tracks of the late optical luminosity lightcurves. Panaitescu & 
Vestrant (2008, 2011) showed some general features of the early optical bumps and plateaus 
in the optical lightcurves. Kann et al. (2010, 2011) compared the optical lightcurves of 
different types of GRBs. 

Different from these previous statistical studies, in a series of papers, we plan to make 
empirical fitting to the observed optical lightcurves and to identify multiple emission com- 
ponents. After decomposing the lightcurves, we plan to perform statistical analyses of the 
parameters of various components, and discuss their physical implications. As the first 
paper in the series, in this paper we present the sample (§2), the general results of the de- 
composition, and a "synthetic" lightcurve that shows eight possible components with distinct 
physical origins (§3). Since flares and the shallow decay component may be directly related 
to the late central engine activity, in this paper we carry out a detailed study of these two 
components: flares in §4 and the shallow decay component in §5, respectively. We discuss 
physical implications in §6, and summarize the results in §7. 
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Data 



We include all the GRBs that have optical afterglow detections before November 2011 
(from Feb. 28 1997 to Nov. 2011) in our sample. A sample of 225 optical lightcurves are 
complied with data reported in the literature. We make an extensive search for the optical 
data from published papers or from GCN Circulars if no published paper is available for some 
GRBs. Well-sampled lightcurves are available for 146 GRBs. In Table 1, we summarize the 
following information for each GRB: redshift, spectral properties of the optical afterglow and 
prompt gamma-ray emission, as well as time intervals of optical observation^^. We collect 
the optical spectral index /3o (defined with the convention Fj, oc z/~'^o)g and the extinction 
Ay of the host galaxy for each burst from the same literature to reduce the uncertainties 
introduced by different authors. Galactic extinction correction is made by using a reddening 
map presented by Schlegel et al. (1998). Since the Ay values are available only for some 
GRBs and the Ay is derived from the spectral fits using different extinction curves, we do 
not make correction for the GRB host galaxy extinction. The ^-correction in magnitude is 
calculated by A; = — 2.5(/3o — 1) log(l + z). For the late epoch data (~ 10^ seconds after the 
GRB trigger), possible flux contribution from the host galaxy is subtracted. The isotropic 
gamma-ray energy {E^^iso) is derived in the rest frame 1 — 10^ keV energy band using the 
spectral parameters. 



3. Lightcurve Fitting and a Synthetic Optical Emission Lightcurve 

The optical lightcurves are usually composed of one or more power-law segments along 
with some flares, humps, or re-brightening features. The mix of different components makes 
the diverse optical afterglow lightcurves. In order to decompose the rich features, we fit 
the lightcurves with a model of multiple components. The basic component of our model is 
either a power-law function 

Fi = Foit-'^ (1) 



or a smooth broken power-law function 



02 



, \ aiuj / -I- \ ^2^' 



ih,lj V^b 



(2) 



^ A full version of the GRB sample with references to the observational data are available in the electronic 
form. 

^An optical spectral index /3o = 0.75 is adopted for those GRBs whose Po is not available. 
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where a, ai, 0:2 are the temporal slopes, tb is the break time, and u measures the sharpness 
of the break. In some afterglow models, a double broken power law lightcurve is expected. 
For example, it is theoretically expected that the afterglow lightcurve may have a shallow 
segment early on due to energy injection, and then transits to a normal decay segment when 
energy injection is over, and finally steepens due to a jet break (e.g. in the canonical X-ray 
afterglow lightcurve, Zhang et al. 2006). We therefore also consider a smooth triple power- 
law function to fit some lightcurves. In this case, we extend equation ([2]) to the following 
function (Liang et al. 2008) 

F^ = {F^'^^ + F^^'')-^''"' (3) 

where UJ2 is the sharpness factor of the jet break at ^{,^2, and 

t ^-"^ 



F, = F,{t^,,) ^—j . (4) 

We developed an IDL code to make best fits with a subroutine called MPFITlj. The 
parameter u is usually fixed to 3 or 1 in our fitting. The approach of our lightcurve fitting is 
as follows. Initially, we introduce a minimum number of components based on eye inspection 
of the global feature of the lightcurve. If the reduced x? is much larger than 1, we continue 
to add more components and re-do the fit, until the reduced Xr becomes close to 1. The 
reduced x^ of & simple power-law fits to the lightcurves of some GRBs, such as GRBs 
020813, 030328, 050416A , and 070110, are ~ 1. However, a smooth broken power-law can 
significantly improve the fito, which reduces > 50% of the Xr value. We therefore adopt 
the smooth broken power-law fit for these GRBs. The Xr values for some lightcurves are 
much smaller than 1, indicating that some model parameters are poorly constrained. For 
these cases, we fix some parameters to make the fits. The erratic fluctuations of some data 
points with small error bars in some GRBs, such as GRB 030329, make Xr much larger 
than 1. We do not add additional components for these lightcurves, so that their x"^ values 
remain much larger than 1. The most challenging problem in our fit is to extract severely 
overlapping flares/bumps from the lightcurves. The slopes of these flares/bumps are usually 
quite uncertain. In our fitting, we first set all the parameters free to get the best fit to the 
global lightcurve, and then adjust the rising slopes to ensure that the fitting curve crosses 
the data point around the peak time of each flare/bump. Finally, we fix the rising slopes and 
perform the best fits again. As an example. Figure [1] shows six lightcurves with the best-fit 
multiple components decomposed. 



■^ http://www.physics.wisc.edu/ craigni/icll/fitting.htmL 

"'Note that the hghtcurves of these GRBs are usually poorly sampled in the initial shallow-decay segment, 
leading to uncertainty to reveal the initial shallow-decay segment. 
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Even though individual optical lightcurves can differ significantly, after synthesizing 
many lightcurves, one can come up with a synthetic afterglow lightcurve, as shown in Figj2j 
In contrast to the five-component canonical X-ray lightcurve (Zhang et al. 2006), the syn- 
thetic optical lightcurve includes eight components that may have distinct physical origins. 
These components are: la. prompt optical flares; lb. an early optical flare from the reverse 
shock; II. early shallow decay segment; III. the standard afterglow component (an onset 
hump followed by a normal decay segment); IV. the post jet break phase; V. optical flares; 
VI. rebrightening humps; VII. late supernova (SN) bumps. The components II-V can find 
their counterparts in the canonical X-ray lightcurve. These components can be distinguished 
based on the parameters of our multi-component fits. For example, flares usually have rapid 
rise and fall. We define a flare if the absolute value of its rising and decaying slopes are 
steeper than 2. If flares occur during the prompt emission phase, they are grouped as la 
(prompt optical flares). A reverse shock flare (lb) is a huge flare (or optical flash) that peaks 
slightly after the end of prompt emission (the decay slope can be somewhat shallower than 
2). All the optical flares (V) afterwards are considered "late" (with respect la and lb), even 
though most of them actually happen in the early afterglow phase. 

We define different components based on theoretical guidance. For example, flares are 
defined as features with both steep rising (ai < —2 with the convention F^ oc t~") and 
decaying (0^2 > 2) slopes. A shallow decay component, on the other hand, is defined as a 
segment whose decay is shallower than what is predicted in the constant energy afterglow 
model. As shown in Table 1, the spectral indices of the optical afterglows are usually 
smaller than 1 within error. This is consistent with the standard external shock synchrotron 
radiation model in the spectral regime of t'm < ^o < ^c (with the electron power-law index 
p > 2), where z/m and z/c are the injection and cooling frequencies, respectively. In this 
spectral regime, one has the closure relation oq = 3/3o/2. A shallow decay component (II) 
is then defined by the condition ao < 3/3o/2 (with the convention F^, oc t~"o) within error. 
Similarly, a post-jet-break decay segment IV is defined as the segment beyond a steepening 
break after the normal decay segment (with ao = 3/3o/2)|3. An afterglow onset feature is 
characterized by a smooth hump peaking at less than 1 hour post trigger, which is followed 
by a normal power-law decay component (HI). A rebrightening hump (VI) is similar to the 
early afterglow onset hump but is much later. It differs from optical flares by much shallower 
rise and decay as well as a much smoother peak. The supernova bump (VII) is a special 



^The radiation physics during the jet-break segment would be the same as the normal decay segment, 
but the break separates two different dynamical evolution regimes of the GRB fireball. We therefore define 
the post-jet-break segment as a new component in our analysis, to echo the definition in the canonical X-ray 
afterglow lightcurve (Zhang et al. 2006). 
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late rebrightening peaking at around 1-2 weeks after GRB trigger, which usually shows a 
red color. 

After decomposing the lightcurves, we are able to group all the identified components to 
one of these eight components. The early optical afterglow lightcurves (t < 10^ s) of about 
one-third GRBs show a smooth hump. Another one-third lightcurves start with a shallow 
decay segment. Twenty-four optical flares are observed in 19 GRBs. Late rebrightenning 
humps are observed in 30 GRBs. A jet like break is detected in 10 GRBs. A clear SN bump 
is detected in 18 GRBs. The detected fraction of each component is marked in the synthetic 
lightcurve (Figj2]). 

We report our statistical results of various components in a series of papers. As the 
first paper of the series, this paper focuses on the optical flares and the shallow decay 
segment. The reason to discuss them together is because they are both likely related to late 
central engine activity of GRBs (see §6 for more discussion). Notice that the prompt optical 
flares and early reversed shock flares are not included in this paper and we'll discuss them 
separately. Throughout the paper, we mark the parameters of the flares and the shallow 
decay segment with the superscripts "F" and "S" , respectively. 



4. Flares 

We get 24 flares in 19 GRBs. A flare is clearly seen in 14 out of the 19 GRBs, as shown 
in Figure [3] along with our best fit results. Some flares may be embedded in the lightcurves 
as shown in Figure [U Notice that most of the well-sampled optical lightcurves are in the 
R band. For a few GRBs, the flares are well-sampled in other bands. We correct these 
lightcurves to the R band with the optical spectral indices. The fitting parameters (the 
flux at peak time and the rising and decaying slopes) of the flares and the derived temporal 
properties, including the peak time (tp), the width {w^) measured at the full- width-half- 
maximum (FWHM), the rising timescale tf , the decay timescale (t^), the ratio of the t^ to 
t^ (-Rfd) ^^d the ratio of the tf to t^ (-Rfp) derived from the fitting parameters are summarized 
in Table 2. With the fitting parameters, we calculate the isotropic peak luminosity {L^ jg^) 
and the total energy release (-Er^sq) i^ ^^^ -^ band. The -Eriso is integrated from t^/b to 
5tp. Our results are reported in Table 2. 

We show the distributions of -Rfp, R^^, t^, w^, L^^^^^ in Figure ID The distribution of 
Rf is clustered around 0.1 — 0.3, being consistent with the expectation of an internal origin 
of the flares. The ratio i?fj is similar to that observed in GRB pulses, but the rising wing of 
some flares are even longer (in log scale) than the decaying wing. The t^ ranges from ~ tens 



of seconds to ~ 10^ seconds. The w^ values are in the same range as t^. The L^^i^o ranges 
from 10^'^ to 10^^ erg s~^, with a typical value of 10^^ erg s~^. 

Relations of w^ and L^ j^^ of the flares as a function of t^ are shown in Figure O A 
tight correlation between w^ and t^ is found. The best fit gives log w^ = —0.32 + 1.01 logtp, 
i.e., w^ ~ '^p/2. The L^^ is anti-correlated with t^ in the burst frame, i.e., logLp^jso4g = 
(1.89 ± 0.52) — (1.15 ± 0.15) log[tp/(l + z)] with a Spearman correlation coefficient 0.85 and 
a chance probability p < 10^^. Therefore, later flares tends to be dimmer and wider than 
earlier flares. 

The possible relations between the flare properties and E^^iso are shown in Figure [61 
Since more than one flares are detected in a few GRBs, we select only the brightest one for 
our analysis. As shown in Figure El -E^^jso is typically smaller than 1/100 of -E^,iso- The flare 
R-band luminosity L^ ^^^ is correlated with gamma-ray luminosity L^^iso, i-e., log L^ jg^/lO"^^ = 
(-3.97±0.60) + (1.14±0.27) logL^,iso/10^° with a Spearman correlation coefficient of r = 0.75 
and a chance probability p ~ 10"^. The flares in GRBs 050401, 060926, and 090726 are 
out of the 3cr region of the fit. Without considering the flares in these three GRBs, it is 
found that the t^ is also tightly anti-correlated with -E^^iso, i-e., logtp^ = (5.38 ± 0.30) — 
(0.78 ± 0.09) logii^^^so/lO^'' (with r = 0.92). Similarly, a tight anti-correlation between Lr^p 
and tp in the burst frame is found without considering the flares in the three GRBs, e.g., 
log[tp/(l -|- z)] = (7.57 ± 0.60) — (1.35 ± 0.17) log £"^,180,50 , with a Spearman correlation 
coefficient 0.91. These results indicate that the optical flares of a GRB that has a larger 
E^^iso tends to peak earlier and brighter. 

It is interesting to study whether optical flares are associated with X-ray flares. Early 
flares are frequently seen in X-ray afterglow lightcurves (Burrows et al. 2005, O'Brien et al. 
2006). However, as shown above, early optical flares are only observed in the lightcurves 
of GRBs 060210, 060926, 090618, and 090726 in our sample. The fraction of GRBs with 
detected early optical flares is much lower than that for X-ray flares. Among the 19 GRBs 
with optical flare detections, 16 had early Swift XRT observations. Their X-ray afterglow 
lightcurves are also shown in Figure 2. Simultaneous observations with XRT during the 
optical flares are available for GRBs 050401, 060206, 060210, 060607A, 060926, 070311, 
071010A, 071031, 080506, 090618, and 100728B. An X-ray flare that may be associated with 
the optical flare is only observed in GRBs 060926, 070311, and 071010A. The optical flares of 
these three GRBs are lagged behind the corresponding X-ray flares. Measuring the lags with 
the peak time of the flares, we get 196 seconds, 7.7 x 10^ seconds, and 2.45 x 10"^ seconds 
for the flares in GRBs 060926, 070311, and 071010A, respectively. The lag is potentially 
proportional to the peak time of the flares with the three flares. 
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5. Shallow decay segment 

A shallow decay segment is defined with a criterion that the initial decay slope of this 
segment is shallower than 3/3o/2 within error. We get a sample of 39 out of 146 GRBs 
that have such a shallow decay segment. Some examples are shown in Figure [3 Figure 
[H] shows the distributions of the decay slopes {af and a|), the break times (t^), and the 
luminosity at the break (Lf^^^) of our sample. Thirty-one shallow decay segments transit 
to a decay slope of 1 ~ 2.5, and 5 have shallow decay segment followed by a sharp drop 
with a decay slope steeper than 2.5. About half of the shallow decay segments look like 
a plateau, with \af^\ < 0.3. The break time ranges from tens of seconds to several days 
after the GRB trigger, with a typical tp of ~ 10^ seconds. The L^^ typically varies from 
10^'^ to 10^^ erg s~^, and even reaches ~ 10'^^ erg s~^ in a few GRBs with an early break. 
The break luminosity L^^ is anti-correlated with t^, as shown in Figure [91 The best fit 
gives logI/|4g = (1.75 ± 0.22) - (0.78 ± 0.08) log[t^/(l + z)], with a Spearman correlation 
coefficient r = 0.86 and a chance probability p < 10^"^. No correlation between E^:^^^ and t^ 
is observed. 

A shallow decay segment is commonly seen in the well-sampled of X-ray afterglow 
lightcurves detected by Swift XRT (e.g., Liang et al. 2007; Evans et al. 2009), except for a 
few GRBs whose XRT lightcurves decay as a single power-law (Liang et al. 2009). It was 
also reported that the X-ray luminosity at the break time is correlated with the break time 
(Dainotti et al. 2010). We over-plot Lb.iso as a function of tb in the burst frame in Figure [91 
One can observe that optical data share a similar relation to the X-ray data. Since the X-ray 
luminosity is measured in the 0.3-10 KeV energy band and the optical luminosity is in the R 
band, the X-ray luminosity lie significantly above the optical luminosity (the uF^, peak is at 
Uc for slow cooling, which may be close to the X-ray band). The observed photon spectral 
indices of the X-ray spectra are ~ 2 (Liang et al. 2007). Therefore, the X-ray energy spectra 
are fiat and the derived Lb,iso — ^b relation in the 1 KeV band is roughly consistent with that 
derived from the entire X-ray band. 

We examine the chromaticity of the shallow decay segments in the X-ray and optical 
bands. The X-ray observations are available for 17 out of the 34 GRBs. We extract the 
underlying afterglow components 11 and 111 (by removing fiares) for the X-ray and opical 
samples, and compare the parameters af, a|, and t^ of the two samples in Figure [TUl It 
is found that the t^ data points are scattered around the equality line, and a tentative 
correlation between the break times of the optical and X-ray lightcurves is observed, with a 
chance probability of the correlation of ~ 0.15. These is no correlation between the decay 
slopes of the X-ray and optical lightcurves. The decay segment prior to the break times in 
the optical bands tends to be steeper than that in the X-ray band, but the post-break slopes 
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are roughly consistent, except for those aj > 2.5 in the optical bands. 

We integrate the total R-band energy release (-Efliso) in the shallow decay phase from 
10 seconds to t^ post trigger. We show E^^^^ as a function of -E^,iso in Figure [TT] A trend 
of rough proportionality between E§^ jg^ and E^^iso is observed. The Spearman correlation 
analysis shows that the chance probability of the correlation between the two quantities is 
~ 6 X 10"'^. Our robust fit yields logE'^^jg^ = 0.40 + 0.47 log -E^iso- The correlation between 
L^iso and -E^^so is much worse. We get logL^^^ = —5.57 + 1.13 log -E^^iso, with a chance 
probability of 0.16, as shown in Figure [TTl 

Physically, there are two types of shallow-decay segment (plateau) , as observed in the X- 
ray band (Liang et al. 2007). The majority of X-ray plateaus are followed by a normal decay 
with decay index typically around -1. These plateaus are likely of an external shock origin, 
with the shallow decay segment caused by continuous energy injection into the blastwave 
(Rees & Meszaros 1998; Dai & Lu 1998; Sari & Meszaros 2000; Zhang & Meszaros 2001a). 
This scenario has been applied to interpret most X-ray plateaus discovered by Swift (Zhang 
et al. 2006; Nousek et al. 2006). A small fraction of plateaus, first found by Troja et al. 
(2007) in GRB 070110 and studied systematically by Liang et al. (2007), are followed by a 
much steeper decay (index steeper than -3), which cannot be interpreted within the external 
shock model. These plateaus are called "internal plateaus" by Liang et al. (2007), since 
they have to be powered by internal dissipation of a late outflow. Looking at our optical 
shallow decay component sample, most lightcurves have a shallow decay component followed 
by a normal decay segment. Nonetheless, we identify two possible internal plateaus in GRBs 
060605 and 080413B, which show superposition of a normal decay segment and a possible 
internal plateau with sharp drop of slope (see Figure [T2li ^l. As shown in Fig. [12], the early 
optical lightcurves of GRBs 060605 and 080413B are a smooth bump and a normal decay 
segment, respectively, which are consistent with the standard afterglow model. Their late 
lightcurve rapidly decays with a slope of a > 2.5, and then flattens to a level consistent with 
the normal decay. This suggests that the plateau is likely internal and is superposed on the 
external component. We also revealed evidence of such a component from the lightcurves of 
GRB 970508 (ending at ~ 10® seconds with a slope 3.0), 050319 (ending at ~ 490 seconds 
and 3.3 x 10^ with slopes 3 and 2.5, respectively). The plateau end times of these GRBs 
range from tens of seconds to several days after the GRB trigger. 



^An issue of defining such cases is that the possibility of the sharp drop is caused by a flare is not ruled 
out. 
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6. Physical implications 

6.1. Prom Prompt Gamma-Ray Pulses to Late Optical Flares: Global 

Evolution of an Erratic Central Engine 

As shown above, flares are an independent component superimposed onto the afterglow 
component. The observed relations between -E^,iso or L^^g^ and t^ indicate that the prompt 
gamma-ray emission and late optical flare emission could have the same physical origin. The 
temporal evolution from the GRB phase to late optical flares may signal global evolution of 
the erratic GRB central engine. 

The most extensively discussed GRB central engine model is a hyper-accreting black 
hole surrounded by a neutrino-dominated accretion disk or torus (NDAFs, e.g., Popham et 
al. 1999; Narayan et al. 2001; Kohri & Mineshige 2002; Di Matteo et al. 2002; Kohri et al. 
2005; Gu et al. 2006; Chen & Beloborodov 2007; Liu et al. 2007, 2008, 2010). The prompt 
gamma-ray phase and the late X-ray and optical flares in a burst are usually well-separated, 
indicating that the central engine may intermittently eject a series of shells during these 
emission episodes. Random collisions of these shells would make internal shocks or magnetic 
turbulent reconnections, which would result in the observed variability (e.g. Kobayashi 
et al. 1997; Zhang & Yan 2011). With gamma-ray data alone, no signiflcant trend of 
width and intensity evolution was found (Fenimore et al. 1995). However, considering both 
prompt gamma-ray pulse and late X-ray and optical flares, we flnd that these episodes are 
correlated and show clear temporal evolution, as shown in FigsJS] and El We also show the 
w — tp relation for single-pulse GRBs observed with CGRO/BATSE and the X-ray flares 
observed with Swift/XKT in comparison with the optical flares. We flnd that they follow 
the same relation (see also Chincarini et al. 2007, 2010; Margutti et al. 2010). The general 
trend is that late flares/pulses tend to be wider and dimmer. This cannot be caused by 
hydrodynamical spreading of the shells ejected at late times, but demand that the central 
engine is ejecting thicker, and dimmer shells at late times (Maxham &; Zhang 2009). This 
may be interpreted as late flares being produced by clumps at larger radii, so that spreading 
during accretion process would increase the accretion time into the black hole (Perna et al. 
2006; Proga & Zhang 2006). 

Liang et al. (2010) discovered a tight correlation between E^^iso and the initial Lorentz 



7, ISO; 



factor of the GRB ejecta, i.e., Fq = 182(E^,iso/10^^erg)°-^^. Replacing -E^,iso with L 
Lii et al. (2011) got Fq = 264(L^^iso/10^^ergs~^)°-^''', and suggested that the correlation 
may be interpreted as a natural consequence of the interplay between neutrino-annihilation 
luminosity and neutrino mass loading from a NDAF. At even lower accretion rate, the 
neutrino annihilation mechanism would be ineflicient to power a jet (Popham et al. 1999; 
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Fan, Zhang & Proga 2005). It is unclear whether luminosity is always proportional to 
the accretion rate rh (it can be maintained so if the jet is magnetically launched), and 
whether the Liso — T relation can be extended to lower luminosities. If one naively extends 
both correlations to lower luminosities, the general trend of decaying L with time [Ljip oc 
^-i.i5±o.i5-j jg consistent with decreasing m with t as expected in several models (e.g. m oc t~^'^ 
of Cannizzo et al. 1990 and rh oc t~^'^^ of Frank et al. 2002). The L — T relation would predict 
that the Lorentz factors of optical flares are below 10. In order to maintain a low Ep for these 
flares, the standard internal shock model would have difficulty (Zhang & Meszaros 2002), 
and one needs to attribute to photosphere emission (Thompson et al. 2007) or magnetic 
dissipations (e.g. Zhang & Yan 2011) in order to account for the observations. 



6.2. The Shallow Decay Segment as Probe of Late Energy Injection 

In the framework of the GRB fireball models, the shallow decay segment followed by 
a normal decay can be interpreted as a blastwave with continuous energy injection. There 
are two types of energy injection, one related to a long-lasting central engine (Dai & Lu 
1998; Zhang & Meszaros 2001a), and another related to distribution of Lorentz factor in the 
promptly ejected outflow (Rees & Meszaros 1998). The two cases have effectively similar 
predictions and cannot be differentiated from the data. The existence of internal plateaus 
observed in some GRBs suggest that at least for some GRBs, indeed a long-lasting central 
engine is at work. This may be related to spindown of the central engine, either a rapidly 
spinning black hole or a rapidly spinning magnetar. Theoretical modeling also suggests 
that the energy of flares can pile up onto the blastwave and make a shallow decay segment 
(Maxham & Zhang 2009). Overlapping optical flares on the shallow decay segment are 
observed in some GRBs, such as 970508 and 000301C. This is good evidence of two emission 
sites: an internal origin of the optical flares and the external shock origin of the shallow 
decay segment. So an observed shallow decay component can be a probe of the central 
engine activity and energy injection into the blastwave. 

For a long-lasting central engine, one may parameterize the central engine luminosity 
history as L = Lot~'^. The external shock closure relation between the decay slope and the 
spectral slope can be written as a = (g — 1) -|- (2 -|- g)/3o/2 in case of Um < i^o < ^c (Zhang 
et al. 2006). With the observed a and f3o, we derive the q values for these GRBs in this 
spectral regime. The typical q value is 0.5, as shown in in Figure [T^ For a black hole - torus 
system, the wind may be driven by neutrino annihilation or the Blandford-Znajek (1977) 
mechanism. For a neutrino-driven wind, the annihilation luminosity can be estimated as 
logLyp = 43.6 -|- 4.891og(m/O.OlM0) -|- 3.4a*, where m is the accretion rate, a* is the spin 
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parameter of the central black hole (W. H. Lei & B. Zhang 2012, in preparation). Assuming 
that the pre-collapse density profile is p oc r^, the mass enclosed within r increases with 
radius as r^ , where r = 3 + r for r > —3 and r ~ for r < —3. Then, the mass fall- 
back rate onto the disk is given as riif oc ^(2^+"^ +3)/{3-t)_ Following Kumar et al. (2008), 
r = —1.8, we get rhf oc t~°'^. If the fall-back mass is accreted into the black hole, we 
obtain logLj^p oc t~^-^. This is inconsistent with q values for most bursts in our sample. 
For the Blandford-Znajek mechanism (1977), the maximum of the power can be estimated 
with Lbz ~ rhc^. If m oc t~°-^ as discussed above, the decay slope is consistent with the q 
values for most GRBs in our sample. Alternatively, if the long-lasting wind is driven by a 
spinning down magnetar, one has L{T) = Loem/i^ + T/TemY, where T^m is the characteristic 
timescale for dipolar spin-down. The predicted q values are or 2, not consistent with the 
typical value q = 0.5. On the other hand, if the shallow decay is caused by injection of flare 
energies into the blastwave (e.g. Maxham & Zhang 2009), the predicted decay slope depends 
on the energetics and temporal distributions of the flares, and can be more flexible. 

One issue to explain the shallow decay segment with energy injection onto the blastwave 
is the chromatic breaks in the optical and X-ray bands (e.g., Panaitescu et al. 2006; Liang 
et al. 2007), as shown in Figures 171 and [TUl This raises the concern regarding whether the 
X-ray and optical emission are from the same emission component (e.g., Liang et al. 2009; 
Racusin et al. 2008). Notice that the decay slopes after the break in the optical and X-ray 
bands are usually consistent with the expectation of external shock models. This implies that 
the radiation from the two energy bands could share a similar origin. Introducing inverse 
Compton scattering in the X-ray band may cause chromatic behavior between the X-ray and 
optical bands (e.g. Panaitescu & Kumar 2000; Zhang & Meszaros 2001b). Alternatively, 
the X-rays may be emission from a long lasting wind. A long-lasting reverse shock may 
introduce further complications (e.g. Uhm et al. 2012). In general, mixing of different 
emission components may be the reason for the complex chromatic behaviors observed in 
different energy bands. 



7. Conclusions 

We have systematically decomposed the optical afterglow lightcurves for 146 GRBs 
before November 2011 that have good quality optical data. By fitting the lightcurves with 
multiple components, we get a synthetic optical lightcurve that includes 8 components with 
distinct physical origins. We plan to study these components in detail in a series of papers, 
and in this paper we focus on the optical flares and the shallow-decay segment. Our results 
can be summarized as the following: 
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• We obtained 24 optical flare events in 19 GRBs. The tp ranges from several tens of 
seconds to several days post the GRB trigger, and it is tightly correlated with the width 
and peak luminosity of the flares, i.e., w^ ~ t^/2 and logL^;gQ4g = (1.89 it 0.52) — 
(1.15 ± 0.15) log[tp/(l + z)], suggesting that flares peaking later tend to be dimmer 
and wider. The parameters t^ and -Eriso are also corrected with -E^,iso, suggesting that 
GRBs with larger E^^iso tend to have optical flares peaking earlier and being brighter. 

• The fraction of GRBs with detected optical flares and the number of flares in a GRB are 
much smaller than the case of X-ray flares. Among the 19 GRBs with detected optical 
flares, 16 have early Swift XRT observations. Only four cases, i.e. GRBs 970508, 
060926, 070311, and 071010A, show detection of associated X-ray flares. The optical 
flares in the three GRBs are lagged behind the corresponding X-ray flares, similar to 
the spectral lag observed in prompt gamma-ray emission, but the time lag is much 
longer what is observed in the prompt phase. 

• We get a sample of 42 shallow decay segments from 39 GRBs. About half of the 
shallow decay segments look like a plateau, with a decay slope af being smaller than 
0.3. Thirty-two out of the 39 shallow decay segments transit to a decay with a slope 
of 1 ~ 2.5, and 5 of them are followed by a sharp drop with a decay slope steeper 
than 2.5. The break times range from tens of seconds to several days after the GRB 
trigger, with a typical t^ ~ 10^ seconds.. No clear correlation between Ef^^^^ and -E^,iso 
is found. 

• The break times of the shallow decay segment in the optical and X-ray bands are 
chromatic for most GRBs, but they are tentative correlated. The decay prior to the 
break time in the X-ray band tends to be steeper than that in the optical band, and 
the decay slopes post the break time in the two energy bands are roughly consistent 
with each other. The L^^go i^ anti-correlated with tf, which is similar to the case for 
X-ray plateaus. 

We discussed the physical implications of the optical flares and the optical shallow decay 
segment, both are related to late GRB central engine activities. The observations strengthen 
the trend that the GRB central engine dies out gradually with the decreasing luminosity with 
time. The late central engine activity can be either erratic (for flares) or steady (for internal 
plateaus), both could add energy to the blastwave to make a shallow decay segment in the 
lightcurve. The observed afterglow is a mix of various emission components of the external 
and internal origins, and the variation of the strengths of different components lead to diverse 
chromatic afterglow behaviors. 
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Table 1. Properties of the GRB sample with well-sampled optical lightcurves 



GRB 


z 


^o 


^v 


rp a 

-t start 


-'end 




^P^ 


a. 


/3 


Refs.'i 


970228 


0.695 


0.78±0.022 


0.5 


69.98 


3180.00 


254.9±23.2 


195±64 


-1.54±0.08 


-2.5±0.4 


(1),(1),(1),(2) 


970508 


0.835 


1.11 




25.57 


7420.00 


42.6±7.1 


145±43 


-1.71±0.1 


-2.2±0.25 


(l),(l),(-),(3) 


971214 


3.42 


0.87±0.13 


0.43±0.08 


46.66 


304.99 


3469.1±354.8 


685±133 


-0.76±0.1 


-2.7±1.1 


(1),(1),(1),(2) 


980326 


1 


0.80±0.4 





36.46 


1850.00 


21.8±4.4 


71±36 


-1.23±0.21 


-2.48±0.31 


(4),(1),(1),(3) 


980425 


0.0085 




1.9±0.1 


564.42 


3880.00 


(83 ±8.7) X lO-'^ 


119±24 


-1±0.3 


-2.1±0.1 


(l),(-),(5),(2) 


980519 




1.07±0.12 




29.98 


231.38 










(-),(6),(-),(-) 


980613 


1.096 


0.6 


0.45 


59.44 


3440.00 


34.8±7.0 


194±89 


-1.43±0.24 


-2.7±0.6 


(1),(1),(1),(3) 


980703 


0.966 


1.013±0.016 


1.5±0.11 


81.26 


343.92 


685.6±59.6 


502±100 


-1.31±0.14 


-2.39±0.26 


(1),(1),(1),(2) 


990123 


1.6 


0.75±0.07 





0.02 


0.61 


62818.1±8375.7 


2030±160 


-0.89±0.08 


-2.45±0.97 


(1),(1),(1),(2) 


990510 


1.619 


0.55 


0.22±0.07 


12.44 


340.24 


1979.4±208.4 


423±42 


-1.23±0.05 


-2.7±0.4 


(1),(1),(7),(2) 


990712 


0.434 


0.99±0.02 


0.5±0.1 


15.25 


320.24 


77.3±3.6 


93±15 


-1.88±0.07 


-2.48±0.56 


(1),(1),(5),(2) 


991208 


0.706 


0.75±0.03 


0.05 


179.52 


613.24 


2230.0 


313±31 






(7),(8),(9),(10^ 
(4),(6),(-),(2) 


991216 


1.02 


0.57±0.08 




38.75 


9470.00 


6534.7±687.9 


641±128 


-1.234±0.13 


-2.18±0.39 


000301C 


2.03 


0.7 


0.09±0.04 


129.17 


4200.00 










(1),(1),(1),(-) 


000418 


1.12 


0.75 


0.96 


214.27 


3930.00 


910.0 


284±21 






(1),(1),(1),(-) 


000630 








73.44 


336.10 




216±56 


-0.67±0.38 


-2.18 


(-),(-),(-),(ll) 


000926 


2.07 


1.00±0.18 


0.18±0.06 


74.48 


505.21 


2710.0±2964.5 


310±20 






(i),(i),(i),(io) 


010222 


1.48 


1.07±0.09 





13.09 


186.21 


8570.9±125.8 


291±43 


-1.05±0.16 


-2.14±0.58 


(1),(1),(1),(11) 


011121 


0.36 


0.8±0.15 





37.35 


1490.00 


780.0 


>700 




>-2 


(1),(1),(1),(12) 


020405 


0.69 


1.43±0.08 





85.04 


882.60 


1328.9±125.7 


364±73 


0.25 


-1.87±0.23 


(1),(1),(1),(12) 


020813 


1.25 


0.85±0.07 


0.14±0.04 


6.05 


5000.00 


7626.5±762.6 


211±42 


-1.05±0.11 


-2.3 


(1),(1),(1),(12) 


020903 


0.251 






57.02 


4548.00 


1.137±0.858 


2.7 


-1 


<-2 


(13),(-),(-),(13) 


021004 


2.335 


0.39 


0.3 


21.12 


2030.00 


499.5±115.3 


267±117 


-1±0.2 




(i),(i),(i),(io) 


021211 


1.01 


0.69 





0.13 


8.96 


111.5±10.1 


47±9 


-0.85±0.09 


-2.37±0.42 


(1),(1),(1),(12) 


030226 


1.98 


0.7±0.03 


0.53 


17.34 


353.70 


1053.4±98.8 


108±22 


-0.95±0.1 


-2.3 


(1),(1),(1),(12) 


030323 


3.37 


0.89±0.04 


0.13±0.09 


34.68 


895.75 


280.0 








(1),(1),(7),(-) 


030328 


1.52 


0.36±0.45 


0.05±0.15 


4.90 


227.47 


2705.8±208.1 


110±22 


-liO.ll 


-2.3 


(1),(14),(14),(12) 


030329 


0.17 


0.5 


0.3±0.03 


11.17 


2860.00 


154.6±14.1 


68±2 


-1.26±0.02 


-2.28±0.05 


(1),(1),(1),(12) 
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GRB 


z 


/3o 


Av 


Tstart^ 




^iso"^ 


^P^ 


a 


/? 


Refs.'i 


030418 








0.29 


7.19 










(-),(-),(-),(-) 


030429 


2.65 


0.75 


0.34 


12.53 


574.04 


216.0 








(1),(1),(1),(-) 


030723 


0.4 


0.66±0.21 


0.32±0.22 


15.00 


6050.00 


4.5 


4.8 


-1 


-2 


(-),(14),(14),(15) 


030725 




2.9±0.6 




335.23 


1280.00 










(-),(16),(-),(-) 


040924 


0.859 


0.7 





0.95 


62.99 


141.3±11.5 


102±35 


-1.17±0.05 




(1),(1),(5),(17) 


041006 


0.716 


0.55 





0.23 


5580.00 


1082.1±282.3 


108±22 


-1.37±0.14 




(1),(1),(1),(17) 


041218 








8.77 


18.62 










(-),(-),(-),(-) 


041219A 


0.31 




6.8±1.6 


0.44 


186.56 










(18),(-),(18),(-) 


050319 


3.24 


0.74±0.42 


0.05±0.09 


0.04 


994.12 


1186.6±187.4 


45±43 


-2±0.2 




(1),(7),(19),(20) 


050401 


2.9 


0.39±0.05 


0.65±0.04 


3.46 


1120.00 


6469.4±1362.0 


119±16 


-0.83±0.13 


-2.37±0.09 


(1),(21),(21),(22) 


050408 


1.2357 


0.28±0.33 


0.73±0.18 


3.35 


3670.00 










(7),(23),(23),(-) 


050416A 


0.65 


1.3 


0.7 


4.09 


32.61 


16.2±1.7 


28.6±8.3 


-1.01 


-3.4 


(5),(24),(5,(17) ^ 


050502A 


3.793 


0.76±0.16 





0.05 


17.85 










(1),(5),(5),(-) ^ 


050525A 


0.606 


0.97±0.1 


0.25±0.16 


0.07 


35.64 


379.5±73.0 


127±5.5 


-l.OliO.ll 




(1),(1),(1),(17) 


050603 


2.821 


0.2±0.1 




34.09 


219.71 


5665.8±323.8 


1333±107 


-0.79±0.06 


-2.15±0.09 


(25),(26),(-),(17) 


050721 




1.16±0.35 




1.48 


248.60 


460.0±90.0 


63±21 


-1.8±0.2 




(-),(27),(-),(20) 


050730 


3.969 


0.52±0.05 


0.12±0.02 


0.07 


72.70 


900.0±300.0 


196±87 


-1.4±0.1 




(1),(5),(21),(20) 


050801 


1.56 


1±0.16 


0.3±0.18 


0.02 


21.65 


46.1±8.4 


44±42 


-1.9±0.2 




(25),(28),(5),(20) 


050820A 


2.612 


0.72±0.03 


0.07±0.01 


0.23 


663.30 


15923.6±1244.0 


1325±277 


-1.12±0.14 




(1),(5),(19),(17) 


050824 


0.83 


0.4±0.04 


0.15±0.03 


0.63 


8990.00 


15.0±4.0 


13±12 




-2.9±0.4 


(21), (21), (21), (20) 


050904 


6.29 


1.31±1.2 


1 


11.05 


459.36 


83951.5±8395.2 


3178±1094 


-1.11±0.06 


-2.2±0.4 


(29),(29),(30),(17) 


050922C 


2.198 


0.51±0.05 





0.74 


606.01 


574.9±195.7 


417±118 


-0.83±0.26 




(31),(5),(5),(17) 


051021 








1.61 


35.82 




99±32 


-0.4±0.8 




(-),(-),(-),(20) 


051028 


3.7 


0.6 


0.7 


8.21 


58.26 


1868.2±168.1 


298±73 


-0.73±0.22 




(-),(32),(33),(33) 


051109A 


2.346 


0.7 




0.04 


13.30 


1198.5±117.5 


539±381 


-1.25±0.5 




(25),(34),(-),(17) 


051111 


1.55 


0.76±0.07 


0.2±0.1 


0.03 


7.59 


1507.3±150.7 


447±175 


-1.22±0.09 


-2.1±0.27 


(25),(53),(-),(22) 


051221A 


0.5459 


0.64±0.05 




11.12 


445.12 


35.5±1.1 


390±190 


-1.34±0.06 




(20),(35),(-),(20) 


060110 


<5 


0.8 




0.03 


4.78 


320.0±60.0 


135±47 


-1.58±0.08 




(25),(36),(-),(20) 
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GRB 


z 


/3o 


Av 


T^ a 

-t start 


-^cnd 


-'^ISO 


^P^ 


a 


/? 


Refs.'i 


060111B 




0.7 


3.6±0.5 


0.03 


13700.00 


1100.0±500.0 


540±280 


-0.9±0.2 






(-),(37),(38),(20) 


060117 








0.13 


0.50 


3600.0±200.0 


72±5 


-1.4±0.1 






(-),(-),(-),(20) 


060121 


4.5 






7.14 


3120.00 


249.6±23.3 


134±32 


0.82±0.38 






(39),(-),(-),(39) 


060124 


2.296 


0.73±0.08 


0.05±0.26 


3.34 


1980.00 


4987.8±733.5 


636±162 


-1.48±0.02 






(25),(40),(7),(17) 


060206 


4.048 


0.73±0.05 


0.01±0.02 


2.89 


201.58 


1568.8±336.2 


381±98 


-1.06±0.34 






(37),(5),(37),(17) 


060210 


3.91 


0.37 


1.18±0.1 


0.06 


7.19 


4150.0±570.0 


575±186 


-1.12±0.26 






(1),(37),(19),(17) 


060218 


0.0331 




0.5 


0.25 


2850.00 


0.5±0.0 


4.9±0.3 


-1.622±0.16 






(25),(-),(5),(17) 


060418 


1.489 


0.78±0.09 


0.12±0.05 


0.08 


7.66 


4859.4±1056.4 


572±114 


-1.5±0.15 






(41),(7),(7),(17) 


060512 


0.4428 


0.68±0.05 




0.11 


5.93 


2.0±0.4 


23±20 








(25),(21),(-),(20) 


060526 


3.21 


0.51±0.32 


0.05±0.11 


0.06 


893.55 


606.4±303.2 


105.2±21.1 


-1.1±0.4 


-2.2±0.4 


(1),(7),(7),(17) 


060605 


3.78 


1.06 





0.07 


6.32 


283.0±45.0 


490±251 


-1±0.44 






(42),(42),(42),(17) 


060607A 


3.082 


0.72±0.27 





0.07 


14.73 


2341.8±148.6 


575±200 


-1.09±0.19 






(37),(7),(19),(ia 


060614 


0.125 


0.47±0.04 


0.11±0.03 


1.55 


1280.00 


21.0±56.2 


55±45 








(21),(21),(21),(ff) 


060714 


2.711 


0.44±0.04 




3.86 


185.46 


1510.0±195.4 


234±109 


-1.77±0.24 






(43),(21),(-),(17) 


060729 


0.54 


0.78±0.03 


0.07±0.02 


0.70 


662.39 


64.9±4.5 


67±25 


-1.8±0.1 






(21),(21),(21),(20) 


060904B 


0.703 


l.lliO.l 


0.08±0.08 


0.02 


163.13 


77.0±9.9 


103±26 


-0.61±0.42 


-1.78±0.23 


(25),(7),(7),(22) 


060906 


3.686 


0.56±0.02 


0.09 


0.66 


13.61 


1726.7±139.3 


209±43 


-1.6±0.31 




(21),(7),(7),(17) 


060908 


2.43 


0.3 


0.05±0.03 


0.83 


7.24 


1334.2±127.9 


124±24 


-0.89±0.2 


-2.24±0.34 


(37),(44),(19),(22) 


060912A 


0.937 


0.62 


0.46±0.23 


1.10 


23.90 


99.8±4.5 


200±110 


-1.7±0.09 




(25),(37),(45),(20) 


060926 


3.2 


0.82±0.01 


0.32±0.02 


0.06 


7.16 


100.0±20.0 


19±18 




-2.5±0.3 


(21),(21),(21),(20) 


060927 


5.6 


0.86±0.03 


0.12 


0.02 


1.17 


5815.2±861.5 


473±116 


-0.93±0.38 




(21),(21),(21),(17) 


061007 


1.261 


0.78±0.02 


0.39±0.01 


0.03 


14.60 


42103.8±4190.4 


498±30 


-0.53±0.05 


-2.61±0.31 


(46),(21),(46),(22) 


061121 


1.314 






7.14 


120.81 


2665±235.3 


606±90 


-1.32±0.05 




(47),(-),(-),(47) 


061126 


1.1588 


0.95 


0.1±0.06 


0.04 


156.38 


28467.5±3272.1 


1337±410 


-1.06±0.07 




(48), (48), (7), (17) 


070110 


2.352 


1.00±0.14 


0.08 


0.66 


34.76 


723±104 


110±50 


-1.57±0.12 




(Swift),(-),(-),(49) 


070125 


1.547 


0.55±0.04 


0.11±0.04 


105.86 


349.05 


8968.2±773.1 


367±51 


-l.liO.l 


-2.1±0.15 


(50),(21),(46),(51) 


070208 


1.165 


0.68 




1.17 


4.85 


28.0±8.0 


60±20 


1 




(25),(37),(-),(20) 


070306 


1.497 


0.7 


5.45±0.61 


12.44 


215.14 


600.0±100.0 


105 


1.67±0.1 






(52),(52),(52),(20) 
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GRB 


z 


iio 


Av 


Tstzxt^ 


-^cnd 


Siso*^ 


^P^ 


a 


/9 


Refs.'i 


070311 




1.0±0.2 


0.8±0.15 


0.07 


350.93 






-1.3±0.1 




(-),(53),(53),(53) 


070318 


0.836 


0.78 


0.44±0.11 


0.06 


87.37 


134.8±32.7 


365±284 


-1.34±0.27 


-2.15±0.36 


(25),(37),(45),(54) 


070411 


2.954 






0.18 


516.63 


1000.0±200.0 


108 


1.7±0.1 




(25),(-),(-),(20) 


070419A 


0.97 


0.8 


0.42±0.37 


0.21 


62.22 


18.7±2.1 


30±7 


0±2 




(25),(24),(7),(20) 


070420 








0.12 


10.84 


3100.0±500.0 


150±40 


1±0.2 




(-),(-),(-),(20) 


070518 


1.16 


0.8 


0.3 


2.11 


311.76 


24.6 


36±33 




-2.1±0.3 


(25),(24),(7),(55) 


070611 


2.04 


0.73 




0.27 


8.87 


44.8±6.6 


67±26 


1 




(25),(37),(-),(55) 


070707 


1.0 


0.75±0.13 




39.52 


3320.00 


8.9±4.2 




-1.19±0.13 




(56),(56),(-),(57) 


071003 


1.605 


1.25±0.09 


0.34±0.11 


0.57 


5.00 


1800.0±600.0 


410±190 


-1.31±0.07 




(58),(58),(7),(55) 


071010A 


0.98 


0.68 


0.64±0.09 


0.32 


523.23 


13.0±2.0 


37±35 




-2.1±0.4 


(25),(59),(7),(55) 


071010B 


0.947 






0.06 


174.46 


173.8±90.0 


52±10 


-1.25±0.49 


-2.65±0.29 


(60),(-),(-),(60) 


071025 


5.2 


0.42±0.08 


1.09±0.2 


0.17 


14.88 


1500.0±300.0 


165±59 


1.67±0.06 




(25),(37),(36),(55) ^ 


071031 


2.692 


0.64±0.01 


0.14±0.13 


0.07 


350.93 


390.0±60.0 


12±11 




-2.3±0.3 


(25),(21),(7),(55) ^ 


071112C 


0.823 


0.63±0.29 


0.23±0.21 


0.13 


69.64 






-1.09±0.07 




(7),(7),(7),(61) 


071122 


1.14 


0.83 




1.30 


9.05 


30.0±10.0 








(25),(37),(-),(-) 


080109 








72.58 


3000.00 










(-),(-),(-),(-) 


080129 


4.349 






0.42 


500.70 










(62),(-),(-),(-) 


080310 


2.4266 


0.42±0.12 


0.19±0.05 


0.30 


124.42 


590.0±100.0 


22±20 




-2.4±0.2 


(25),(7),(7),(55) 


080319A 


<2.2 


0.77±0.02 




0.15 


4.46 


800.0±100.0 


105±35 


-1.6±0.1 




(25),(21),(-),(55) 


080319B 


0.937 






0.01 


4590.00 


52639.3±4024.7 


>1382 


1.09±0.02 




(55),(-),(-),(55) 


080319C 


1.949 


0.77±0.02 


0.59±0.12 


0.08 


1.43 


5206.3±1041.3 


307±58 


-1.01±0.08 


-1.87±0.39 


(21),(44),(7),(22) 


080330 


1.51 


0.49 


0.19±0.08 


0.09 


116.56 


41.0±6.0 


20±19 




-2.4±0.5 


(25),(59),(7),(55) 


080413A 


2.433 


0.67 


0.13±0.07 


1.17 


18.34 


1855.0±397.0 


126±42 


-1.15±0.29 


-2.12±0.33 


(25),(59),(7),(54) 


080413B 


1.1 


0.25±0.07 




0.08 


5190.00 


175.7±21.9 


67±8 


-1.24±0.26 


-2.77±0.22 


(25),(59),(-),(54) 


080506 




0.95±0.05 




0.21 


5.37 


190.0±40.0 


67±28 


-1.70.2 




(-),(63),(-),(55) 


080603A 


1.67842 






0.11 


350.44 










(Swift), (-),(-),(-) 


080710 


0.845 


0.8±0.09 


0.11±0.04 


0.42 


266.59 


80.0±40.0 


300±200 


-1.3±0.2 




(25),(7),(7),(55) 


080721 


2.602 


0.68±0.02 


0.6 


38.40 


559.52 


18915.6±556.3 


485±37 


-0.93±0.05 


-2.43±0.26 


(30),(21),(64),(22) 
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GRB 


z 


/3o 


^v 


Tlstart^ 


^ end 


Siso*^ 




a 


ii 


Rcfs.d 


080804 


2.2 


0.43 




1.16 


26.11 


1600.0±700.0 


410±200 


-1±0.1 




(25),(59),(-),(55) 


080810 


3.35 


0.44 


0.16±0.02 


0.04 


7.90 


3000.0±2000.0 


313.5 


-0.91 




(25),(7),(7),(65) 


080913 


6.7 


0.79±0.03 


-0.58±0.67 


0.58 


870.04 


710.7±88.8 


135±47 


-0.4±0.9 




(25),(21),(7),(55) 


080928 


1.692 


1.08±0.02 


0.29±0.03 


0.39 


13.43 


280.0±50.0 




-1.8 




(21),(21),(21),(65) 


081008 


1.967 






0.11 


184.52 


630±350 


117±50 


-1.26±0.24 




(66),(-),(-),(66) 


081028 


3.038 






11.19 


57.50 


1081.4±1405.9 


67±13 


1.3±0.4 




(25),(-),(-),(55) 


081029 


3.85 


libO.Ol 


0.03±0.02 


0.53 


252.67 


1513.6±558.6 








(67),(67),(67),(-) 


081109A 


0.98 






0.17 


66.60 


530.0±80.0 


99±40 


-1.27±0.34 


-2.19±0.42 


(25),(-),(-),(54) 


081126 








0.10 


0.54 


900.0±200.0 








(-),(-),(-),(55) 


081203A 


2.1 


0.596 


0.09±0.04 


0.08 


5.76 


1700.0±400.0 


201±75 


-1.44±0.06 




(25),(7),(7),(55) 


090102 


1.547 


0.74 


0.12±0.11 


0.04 


264.55 


1400.0±500.0 


370±220 


-1.36±0.1 




(25),(7),(7),(55) 


090313 


3.375 


0.71 


0.34±0.15 


0.20 


7870.00 


460.0±50.0 


55±51 


-1.9±0.3 




(25),(59),(7),(55) ^ 


090323 


3.57 


0.65±0.13 


0.14±0.04 


162.43 


768.96 


33626.4 


697 


-0.89 




(68),(68),(68),(65) T' 


090328 


0.736 


1.19±0.21 


0.22±0.12 


57.89 


1070.00 


1902.6 


653 


-0.93 


-2.2 


(68), (68), (68), (65) 


090426 


2.609 


0.76±0.14 




0.09 


10.75 


42.04.0 


45±43 




-2±0.3 


(69),(70),(-),(55) 


090510 


0.903 






0.11 


103.79 




3900±280 


-0.58±0.06 


-2.83±0.20 


(71),(-),(-),(72) 


090618 


0.54 


0.5 


0.3±0.1 


0.08 


72.58 


2476.6 


155.5 


-1.26 


-2.5 


(73),(73),(73),(65) 


090726 


2.71 






0.20 


3.02 


186.9±17.2 


27±22 


-1.2±1.3 




(25),(-),(-),(55) 


090812 


2.452 


0.36 




0.03 


0.14 


4585.8±597.4 


190±65 


-1.5±0.3 




(25),(59),(-),(55) 


090902B 


1.822 


0.68±0.11 


0.2±0.06 


4.80 


563.79 










(74),(74),(74),(-) 


090926A 


2.1062 


0.72±0.17 


0.13±0.06 


73.16 


102.15 


26562.4±963.4 


412±20 


-0.74±0.01 


-2.34±0.01 


(7),(7),(7),(75) 


091029 


2.752 


0.57 




14.26 


43.09 


849.5±35.4 


61.4±17.5 


-1.46±0.27 




(25),(59),(-),(76) 


091127 


0.49 


0.43±0.10 





2 


7.93 


8820.00 


152.6±7.5 


21.3±3 


-1.95±0.1 




(77), (77), (77), (78) 


100219A 


4.6667 








0.94 


35.0 


359±64.3 


140±0 


-1.34±0 




(79),(-),(-),(79) 


100316D 


0.059 








40.26 


4960.00 










(80),(-),(-),(-) 


100418A 


0.6235 








1.01 


1371.57 










(Swift),(-),(-),(-) 


100728B 










0.16 


5.64 






1.55±0.14 




(-),(-),(-),(81) 


100901A 


1.408 








0.64 


543.01 


245.50 




1.52±0.21 




(82),(-),(-),(82) 



Table 1 — Continued 



GRB 



/3o A^ 



T., 



T, 



end 



Eu 



E^ 



/? 



Rcfs; 



100906A 


1.727 ... 


0.05 


10.94 






. (Swift), (-),(-), (- 


101024A 




0.22 


160.44 


56.25±5.54 


-1.4±0.8 


■ (-),(-),(-),(83) 


101225A 


0.33 


5.72 


3499 




-1.8±0.32 . 


. (84),(-),(-),(84) 


110205A 




14.34 


384.19 


230±65 


-0.59±0.06 . 


■ (-),(-),(-),(85) 


110213A 


1.46 


0.10 


183.37 720.0±10.0 


98.4±6.9 


-0.44±0.05 . 


. (85),(-),(-),(85) 


110918A 




... 122.43 


1410 






• (-),(-),(-),(-) 
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(12) Ghirlanda et al.(2004); (13) Sakamoto et al.(2004); (14) Kann et al.(2006); (15) Butler et al.(2005); (16) Pugliese et al.(2005); (17) Nav£^ et 
al.(2008); (18) Gotz, D. et al (2011) (19) de Ugarte Postigo et al.(2011); (20) Butler et al.(2007); (21) Zafar et al.(2011); (22) Ukwatta et al.(2(g)); 
(23) de Ugarte Postigo et al.(2007); (24) Xin et al.(2010); (25) Robertson, Brant et al. (2011); (26) Grupe et al.(2006); (27) AntonelH et al.(2006); 
(28) de Pasquale et al.(2007); (29) Kann ct al.(2007); (30) Berger et al.(2007); (31) Price et al.(2006); (32) Urata et al.(2007); (33) Castro-Tirado 
(2006); (34) Yost et al.(2007); (35) Soderberg et al.(2006) (36) Perley et al.(2009); (37) Fynbo et al.(2009); (38) Klotz et al.(2006); (39) Donaghy 
et al.(2006); (40) Misra et al.(2007); (41) Prochaska et al.(2007); (42) Ferrero et al.(2009) (43)Jakobsson et al.(2006); (44) Covino et al.(2010); (45) 
Schady et al.(2011); (46) Schady et al.(2008); (47) Golenetskii et al(2006); (48) Perley et al.(2008); (49) Troja et al.(2007); (50) De Cia et al.(2011); 
(51) Bellm et al.(2008); (52) Jaunsen et al.(2008); (53) Guidorzi et al.(2007); (54) Krimm et al.(2009); (55) Butler et al.(2010); (56) Piranomonte et 
al.(2008) (57) McGlynn et al.(2008); (58) Kriihler et al.(2009); (59) Greiner et al.(2011); (60) Golenetskii et al. (2007); (61) Krimm et al. (2007); (62) 
Greiner et al (2009) (63) Uehara et al.(2010); (64) Starling et al.(2009); (65) Guetta et al.(2011); (66) Yuan et al.(2010); (67) Nardini et al.(2011); 
(68) McBreen (2010); (69) Thonc et al.(2011); (70) Nicuesa et al.(2011); (71) Giuliani et al.(2010); (72) Asano, Katsuaki et al.(2010) ; (73) Cano 
et al.(2011); (74) Pandey, S. B. et al.(2010) (75) Zhang et al.(2011); (76) Barthelmy et al. (2009); (77) Vergani et al.(2011); (78) Golenetskii et al. 
(2009); (79) Mao et al.(2011); (80) Bufano et al.(2011); (81) Barthelmy et al. (2010); (82) Sakamoto et al. (2010); (83) McBreen et al. (2010b); (84) 
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Table 2. Properties of the Flares in Our Sample 



GRB(Band) 


P a 


ai 


a2 


tp 


Lv,K 


EkAso'^ 


w'- 


tr^ 


id^ 


^rd 


-Rrp 


970508(R) 


13G.0±9.7 


-5.00±1.10 


2.50±0.11 


15G.00±10.98 


G.47±G.03 


8.38 


134.73 


45.77 


88.96 


0.51 


0.28 


000301C(R) 


90.8±2.3 


-6.00 


2.11 


138.G6±G.80 


1.95±G.05 


20.18 


85.00 


22.37 


62.63 


0.36 


0.16 


000301C(R) 


43.1±2.4 


-3.20 


11.56 


336.76±2.12 


G.93±G.05 


37.93 


101.27 


69.04 


32.23 


2.14 


0.21 


021004(R) 


125.4±14.4 


-3.00 


3.00 


68.90±2.5G 


2.54±G.29 


30.67 


43.07 


18.25 


24.82 


0.74 


0.26 


050401 (R) 


61.1±14.0 


-5.48 


8.94 


24.58±0.82 


1.90±G.44 


5.10 


6.42 


3.50 


2.92 


1.20 


0.14 


051109A(R) 


623.8±75.9 


-2.10 


2.24±0.56 


7.61±G.59 


18.52±2.25 


20.54 


6.60 


2.67 


3.93 


0.68 


0.35 


060121(R) 


7.0±0.8 


-5.00±2.78 


4.GG±0.21 


34.GG±0.46 


O.G6±O.GG7 


0.65 


14.53 


5.99 


8.53 


0.70 


0.17 


060206(R) 


2246.G±39.8 


-9.74 


2.21±0.11 


3.45±G.01 


229.56±4.07 


25.37 


1.75 


0.35 


1.40 


0.25 


0.10 


060206(R) 


593.7±34.5 


-5.70 


2.5G 


6.65±0.05 


60.68±3.52 


20.54 


3.45 


1.02 


2.43 


0.42 


0.15 


060210(R) 


888.9±245.3 


-3.20 


5.72 


0.66±G.03 


47.54±13.12 


4.75 


0.28 


0.15 


0.13 


1.17 


0.23 


060607A(H) 


856.2 


-6.08±2.67 


2.41±0.28 


4.26±G.39 


47.66 


22.81 


2.49 


0.73 


1.76 


0.42 


0.17 


060607A(H) 


1061.7 


-3.00 


10.89±2.77 


2.17±G.04 


59.10 


24.85 


0.72 


0.48 


0.24 


1.96 


0.23 to 

CO 
0.28 


060926(V) 


1115.3±215.9 


-2.55±1.08 


3.49±0.77 


G.63±G.05 


77.80±15.06 


3.69 


0.36 


0.17 


0.19 


0.93 


060926(V) 


2629.7±481.9 


-3.47±G.87 


2.GG 


G.09±G.01 


183.45±33.62 


1.06 


0.07 


0.02 


0.05 


0.48 


0.24 


07031 1(R) 


260.8±42.4 


-5.10±1.19 


5.31±1.2G 


2G6.89±9.40 






90.50 


40.72 


49.78 


0.82 


0.20 


071010A(R) 


157.9±5.4 


-2.08 


2.39 


80.G7±1.4G 


0.63±0.02 


10.84 


66.30 


27.80 


38.50 


0.72 


0.35 


071025(J) 


1650.9±119.1 


-5.67±1.67 


3.00 


1.71±G.04 






0.86 


0.29 


0.56 


0.52 


0.17 


071031(R) 


16.8±0.1 


-8.27 


3.00 


16.57±0.G2 


0.63±0.01 


0.58 


7.35 


2.17 


5.19 


0.42 


0.13 


071031(R) 


49.8±0.3 


-3.51 


3.40 


6.36±G.01 


1.86±0.01 


1.30 


3.45 


1.47 


1.98 


0.74 


0.23 


080506(R) 


405.6±98.4 


-6.84±7.49 


2.38±2.49 


1.19±G.13 






0.68 


0.19 


0.49 


0.39 


0.15 


090313(R) 


1452.5±190.5 


-6.43±1.24 


2.GG 


19.83±0.69 


97.14±12.74 


107.35 


13.39 


3.43 


9.96 


0.34 


0.17 


090618(R) 


33652.9±1055.6 


-5.00 


2.23±0.G7 


G.14±G.OO 


31.34±0.98 


0.68 


0.09 


0.03 


0.06 


0.44 


0.19 


090726(R) 


873.1 


-2.21±G.69 


5.GG 


0.52 


38.08 


2.49 


0.28 


0.15 


0.12 


1.27 


0.31 


100728B(R) 


64.5±16.1 


-8.00 


5.6G±2.29 


3.00±G.39 






1.71 


0.70 


1.01 


0.69 


0.22 



'^In units of 10^^^ erg cm^^ s^^. 
''In units of kilo seconds. 
■^In units of 10^'"' erg. 
''in units of 10'** erg. 



00 

o 



Table 3. Properties of the Shallow Decay Segments in Our Sample 



GRB 


Fm" 


ai 


^2 


f ^ 


ip,R' 


^Raso"^ 




q 


aiiPo) 


970508(R) 


23.8±3.9 


0.10±0.04 


3.00±0.09 


1000.00±58.24 


85±14 


70.0 


114.3 


-0.006 


1.665 


000301C(R) 


19.0±1.1 


0.20±0 


3.42 


602.81±16.55 


408±24 


172.9 


365.8 


0.370 


1.05 


010222(R) 


323.7±39.7 


0.47±0.06 


1.25±0.07 


32.07±4.89 


4754±583 


157.6 


313.2 


0.261 


1.605 


020813(R) 


348.9±177.3 


0.50±0.35 


1.37±0.11 


29.91±14.82 


2230±1130 


90.5 


141.8 


0.453 


1.275 


021004(R) 


1841.4±88.9 


0.28±0.02 


1.50±0.05 


9.91±0.56 


37249±1798 


564.0 


589.7 


0.742 


0.585 


030226(R) 


54.0±4.7 


0.70±0.04 


2.92 


102.35±5.17 


1097±95 


142.4 


343.0 


0.741 


1.05 


030328(R) 


245.4±66.1 


0.41±0.13 


1.28±0.06 


16.11±4.74 


1999±539 


65.7 


60.9 


0.893 


0.54 


030429(R) 


18.4±1.6 


0.86±0.03 


3.53±0.00 


218.47±8.85 


606±53 


111.8 


582.7 


0.805 


1.125 


030723(R) 


41.4±1.6 





2.09 


103.16±2.52 


21±1 


3.5 


2.6 


0.256 


0.99 


040924(R) 


833.4±127.0 


0.05±0.32 


1.29±0.03 


1.72±0.32 


2484±379 


6.9 


5.4 


0.257 


1.05 


041006(R) 


325.4±16.9 


0.42 


1.27±0.01 


11.35±0.53 


586±30 


13.6 


12.6 


0.683 


0.825 


050319(R) 


55.0±4.8 


0.40±0.03 


1.40±0.08 


121.11±15.25 


3516±307 


381.3 


798.3 


0.479 


1.11 


050319(R) 


680.2±54.7 


0.32±0.05 


2.07 


0.54±0.04 


43476±3495 


16.7 


35.9 


0.424 


1.11 


050408(R) 


19.8±0.8 


0.52 


1.39 


40.39±2.62 


98±4 


8.6 


8.8 


1.091 


0.42 


050416A(R) 


27.0±1.6 


0.39 


1.32 


15.33±1.11 


57±3 


1.2 


1.6 


0.053 


1.95 


050730(R) 


590.0±83.0 


0.33±0.04 


1.73±0.06 


11.36±1.73 


41271±5806 


549.0 


795.1 


0.641 


0.78 


050801 (R) 


8175.1±1158.8 


0.10±0.13 


2.41±0.30 


0.29±0.01 


128276±18183 


25.8 


45.4 


0.064 


1.5 


050801 (R) 


1647.6±292.7 


0.20±0.28 


1.36±0.05 


1.51±0.30 


25852±4593 


40.5 


56.0 


0.132 


1.5 


050922C(R) 


1912.5±385.0 


0.18±0.07 


1.58±0.04 


4.08±0.91 


39381±7928 


260.4 


230.6 


0.535 


0.765 


051021(R) 


356.5±329.6 


0.26±0.68 


1.50 


4.93±3.20 








0.371 


1.125 


051109A(R) 


4945.0±2565.0 


0.40±0.08 


1.04±0.04 


0.44±0.32 


146835±76164 


89.3 


107.0 


0.518 


1.05 


051111(R) 


2560.0±469.0 


0.81±0.01 


2.1±0.7 


3.06±0.66 


25G7G±4590 


96.1 


244.6 


0.761 


1.14 


060210(R) 


623.8±127.8 


0.07±0.09 


1.21±0.05 


0.96±0.20 


33362±6834 


46.2 


40.4 


0.591 


0.555 


060526(R) 


334.9±19.7 


0.58±0.03 


1.82±0.03 


24.95±1.18 


15070±886 


441.7 


903.8 


0.855 


0.765 


060605(R) 


299.2±51.4 


0.18±0.12 


3.74±0.56 


23.28±2.20 


44198±7597 


299.1 


1486.9 


0.077 


1.59 


060729(U) 


4441.9±1119.1 


0.13±0.17 


2.65±1.66 


4.13±0.97 


4669±1176 


14.3 


26.4 


0.253 


1.17 


060927(V) 


21362.0±11727 


0.38±0.63 


2.63±0.97 


0.05 


(6.8 ±3.6) X 10^ 


99.4 


338.6 


0.287 


1.29 


061126(R) 


179.4±36.4 


0.40±0.12 


1.29±0.05 


14.81±3.03 


1299±263 


25.2 


35.9 


0.308 


1.425 



00 









Table 3 — Continued 














GRB 


P a 


ai 


a2 




ip,R' 


EK,iso'^ 




q 


MM 




070110(U) 


228.1±162.3 


0.16±G.27 


0.98±G.78 


14.29±19.66 


2302±1638 


4G.6 


46.5 


0.477 


1.5 




07041 1(R) 


32.2±1.1 


G.49±G.01 


1.90 


108.99±2.26 


17G2±57 


180.5 


398.5 


G.541 


1.125 




070518(R) 


9.8±G.7 


G.65 


1.90 


30.00 


64±4 


3.G 


5.2 


G.607 


1.2 




070707(R) 


4.3±G.9 


0.44±G.14 


3.20±G.24 


12G.01±13.81 








0.502 


1.125 




071003(R) 34500.0±3520.0 


0.89±G.04 


1.86±G.02 


G.18±G.01 


585510 


5974G 


88.4 


G.392 


1.875 




071010A(R) 


337.3±59.5 


G.29±G.17 


1.55 


12.53±3.1G 


1347±238 


25.6 


27.6 


G.459 


1.02 




080413A(R) 31509.9±4714.6 


G.56±G.08 


5.02±1.59 


G.07±G.01 


977554±146264 


54.1 


97.6 


G.670 


1.005 




080413B(R) 


88.4±5.6 


G.04±G.03 


2.04±G.09 


159.G2±9.32 


335±21 


91.0 


67.6 


G.704 


0.375 




081029(R) 


536.4±39.7 


G.71±G.06 


10.38±0.GG 


2.73±G.10 


4747G±3510 


69.0 


107.3 


1.710 


1.5 




090426(R) 


155G.0±132.0 


G.55±G.06 


1.75±G.OO 


G.36±G.03 


4992G±4250 


16.8 


33.6 


G.571 


1.14 




090426(R) 


69.4±7.83 


G.27±G.06 


2.29±G.OO 


30.48±1.69 


224G±252 


5G.2 


108.4 


G.367 


1.14 




090618(R) 


594.0±25.0 


G.62±G.OO 


1.63±G.06 


31.12±1.46 


553±23 


37.9 


44.8 


G.899 


G.75 


00 


091127(1) 


9G2.4±13.6 


G.45 


1.48 


27.47±0.57 


595±9 


41.1 


32.4 


1.161 


0.645 


to 


101225A(R) 


6.9±G.5 


G.15 


1.30 


368.78±33.90 








0.295 


1.125 




^In units of 10 


^''erg cm ""s . 




















''In units of ks. 






















''In units of 10' 


*2 erg s-\ 
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Fig. 1. — Examples of our model fits (solid lines) to the optical lightcurves with multiple 
components (dashed or dash-dotted lines). The sohd lines represent the best fit to the 
data. Simultaneous X-ray data observed with Swift/XBT (crosses with error bars) are also 
presented. 
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Fig. 2. — A synthetic cartoon lightcurve of multiple optical emission components based on 
our analysis. 
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Fig. 3. — Optical afterglow lightcurves with clear detections of at least one optical flare. The 
line styles and symbols are the same as Figure 1. 
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Fig. 4. — Optical flare parameter distributions in our sample. 
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Fig. 5. — Pair correlations between flare parameters in our sample. Lines are the regression 
lines. 
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Fig. 6. — -^Riso ^^^ -^Riso ^s ^ function of -E^,iso for the flares in our sample. Lines in the 
right panel are the regression line and 3 a region. 
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Fig. 7. — Optical afterglow lightcurves with clear detections of the shallow decay segment (s) 
The symbols are the same as Figure 1. 
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Fig. [7] — Continued 
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Fig. [7] — Continued 
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Fig. 8. — Parameter Distributions of the shallow decay segments for the GRBs in our sample. 
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Fig. 9. — L^ iso (left) and R^ jg^ as a function of tf/{l + z) for the GRBs with a shallow decay 
segment in their optical lightcurves. The grey circles are for the X-ray data from Dainotti 
et al. (2010). Lines are the best fit lines. 
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Fig. 10. — Comparisons of the decay slopes and the break times in the optical and X-ray 
bands. The dashed hnes are the equality hnes. 
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Fig. 11. — -E'7,iso as a function of -E'R,iso, -Z^r.iso, and tp/(l + z) for the GRBs with a shallow 
decay segment in their optical lightcurves. The lines are the robust fits to the data. 
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Fig. 12. — Optical afterglow lightcurves with possible detections of an internal plateau. The 
symbols are the same as Figure 1. 
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Fig. 13. — Distribution of the q parameter of our sample. 



